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Abstract: 

With up to 7 million spins, the existence of spontaneous magnetization 
of Ising spins on directed Barabasi-Albert networks is investigated by Monte 
Carlo simulations. We confirm our earlier result that the magnetization for 
different temperatures T decays after a characteristic time t(T), which we 
extrapolate to diverge at zero temperature by a modified Arrhenius law, or 
perhaps a power law. 
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Introduction: 

The Ising magnet is since decades a standard tool of computational 
physics [1]. We apply it here to scale-free networks [2], where previous sim- 
ulations [3] indicated a Curie temperature increasing logarithmically with 
increasing system size N. In contrast to that work we use here directed [4] 
as opposed to undirected networks and then apply the standard Glauber 
kinetic Ising model [1] to the fixed network. We try to improve our previ- 
ous note [5] by making the system an order of magnitude larger, correcting a 
programming error, and comparing networks with two and seven neighbours. 

Directed scale-free network: 

Putting Ising spins onto the sites (vertices, nodes) of a network, we sim- 
ulate our Ising magnetic model on directed Barabasi-Albert networks. The 
Barabasi-Albert network is grown such that the probability of a new site to 
be connected to one of the already existing sites is proportional to the num- 
ber of previous connections to this already existing site: The rich get richer. 
In this way each new site selects exactly m old sites as neighbours. 
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Decay time tau for 7 million with m=2 (left +) and 2 million with m=7 (right +); and smaller systems 
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Figure 1: Characteristic time for M(r) =3/4 using 7 million spins for m = 2 
neighbours and 7 million spins for m = 7 neighbours (+). Ten, hundred, and 
thousand times smaller systems are denoted by x, stars, and squares. We 
plot the median over nine samples in this log- log plot. The two straight lines 
have negative slopes 8 (left) and 12 (right). 

Then each spin is influenced by the fixed number m of neighbours which 
it had selected when joining the network. It is not influenced by other spins 
which selected it as neighbour after it joined the network. 

At each step, a new spin is added which builds m new connections neighb, 
randomly to already existing spins. The probability for an existing spin to 
be chosen as neighbour is proportional to the number of its neighbours. 

Ising simulations: 

The Ising interaction energy is 

e = -jJ2J2 s ^ (s t = ±i) 

i k 

where the sum over k goes only over the m neighbours which site % had 
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Arrhenius tests: 0.09T + 0.03*T*T*T (m=7, right) and 0.7*t + 2*T*T*T (m=2, left) 
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Figure 2: Same data (except for short times) plotted as l/ln(r) versus T. 



selected when it joined the network. We measure the temperature in units 
of the usual Curie temperature of the square-lattice Ising model. 

First we initialize a directed Barabasi- Albert network of N sites with m 
neighbours (all m initial spins are connected with each other and themselves), 
here m = 7. We put an Ising spins onto every site, with all spins up, because 
we test here for ferromagnetism. Then with the standard Glauber (heat bath) 
Monte Carlo algorithm spins we search for thermal equilibrium at positive 
temperature. Time t is measured in Monte Carlo steps MCS per spin. 

As shown earlier [5] the magnetization M(t) as a function of time t shows 
strong fluctuations and sometimes very rapid changes in a very short interval. 
These problems became worse when we corrected a programming error. Now 
the magnetization reduces to a temperature-dependent metastable value after 
the first time step, stays there for a long time, and then flips or tries to 
flip. We follow [5] and define a characteristic time r as that time where the 
magnetization has decreased to 3/4 of its initial value, Fig.l. This double- 
logarithmic plot first suggests a power-law divergence r(T — > 0) — > oo. But 
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we see slight curvature, unusually high exponents near 10, and an exponent 
varying with the number m = 2 or 7 of neighbours. More plausible thus is 
the modified Arrhenius plot of Fig. 2 which suggests r cx exp(const/T) for 
low temperatures, or l/ln(r) oc T + . . . This const in the exponential is about 
1.4 for m = 2 and about 11 for m — 7. 

Conclusion: 

In this way we confirmed the asymptotic Arrhenius extrapolation 1/ In r oc 
T of [5], meaning that at all finite temperatures the magnetization eventually 
vanishes: No ferromagnetism. 
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